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Recent exoplanet surveys have predicted a very large population of planetary systems in our 
galaxy, more than one planet per star on the average, perhaps totalling about two hundred billion. 
These surveys, based on electro-magnetic observations, are limited to a very small neighbourhood of 
the solar system and the estimations rely on the observations of only a few thousand planets. On the 
other hand, orbital motions of planets around stars are expected to emit gravitational waves (GW), 
which could provide information about the planets not accessible to electro-magnetic astronomy. 

The cumulative effect of the planets, with periods ranging from few hours to several years, is 
expected to create a stochastic GW background (SGWB). We compute the characteristic GW strain 
of this background based on the observed distribution of planet parameters. We also show that the 
integrated extragalactic background is comparable or less than the galactic background at different 
frequencies. Our estimate shows that the net background is significantly below the sensitivities of 
the proposed GW experiments in different frequency bands. However, we notice that the peak of the 
spectrum, at around 10“® Hz, is not too far below the proposed space based GW missions. A future 
space based mission may be able to observe or tightly constrain this signal, which will possibly be 
the only way to probe the galactic population of exoplanets as a whole. 

PACS numbers: 


I. INTRODUCTION 

Einstein’s theory of General Relativity predicts that 
massive bodies with oscillating mass-quadrupole moment 
emit Gravitational Waves (GW) [1, 2]. Planets orbiting 
around stars are therefore expected to emit GW. While 
the signal from a single planet may be too weak to detect 
by the current and upcoming GW detectors [3, 4], it was 
important to check if the huge number of planets pre¬ 
dicted by the recent exoplanet missions, like Kepler [5], 
can constitute a detectable stochastic GW background 
(SGWB). SGWB is an incoherent superposition of signals 
from a large set of sources where an individual source can 
not be resolved [6]. SGWBs of different cosmological and 
astrophysical origin have been estimated in literature [6- 
21] and several upper limits have been placed using data 
from various experiments and many sophisticated anal¬ 
ysis techniques [6, 9, 22-32]. In this paper we estimate 
the background arising from the predicted population of 
exoplanets in the Milky Way galaxy due to their orbital 
motion in their planetary system. The current exoplanet 
missions have observed few thousand planets in a small 
neighbourhood of the solar system. Based on these mea¬ 
surements a population model was estimated, predicting 
more than one planet per star, totalling about two hun¬ 
dred billion planets in our galaxy [33]. Using the ob¬ 
served distributions of planet parameters (mass, elliptic- 
ity, orbital period and semi-major axis) from the publicly 


“Electronic address: ainz®iucaa.ernet.in 
fElectronic address: shilpakastha@imsc.res.in 
^Electronic address: sanjit@iucaa.ernet.in 

S Current Affiliation: The Institute of Mathematical Sciences, 
Chennai 600113, India 


available exoplanet databases we estimate the frequency 
spectrum of the SGWB created by all the planetary sys¬ 
tems in our galaxy. We also compute the strength of 
the background from the Andromaeda galaxy and check 
if the extragalactic component, integrated over the rest 
of the universe, has any significant contribution as com¬ 
pared to its galactic counterpart. 

We then explore the possibility of detection of this 
background. Though direct detection of GW has not 
been possible yet, observation of decay in orbital period 
of binary pulsar PSR B1913-I-16 over three decades pro¬ 
vides a convincing evidence of the existence of GW [34, 
35]. Direct detection of GW is one of the most impor¬ 
tant challenges in the current research in Astrophysics 
and Gosmology, which is likely to happen in the next few 
years, perhaps through different windows of observation. 
While the ground based laser interferometric GW obser¬ 
vatories [36-38] are most promising for detecting GW 
sources at ~ lOOHz, a low frequency stochastic back¬ 
ground may be detected soon though the measurements 
of the Gosmic Microwave Background “B-mode” polari¬ 
sation anisotropy [39-41] and the Pulsar Timing Arrays 
(PTA) [42-45] . The proposed space based detectors [46- 
48], planned to be launched in the next decades, are 
expected to observe different GW signals at very high 
Signal-to-Noise Ratio (SNR) compared to the other de¬ 
tectors. For this work PTA and space based detectors 
are the relevant ones, as their frequency bands have in¬ 
tersections with the frequency range of the planet orbits. 
We discuss the possibility of detection of the SGWB from 
exoplanets using these two kinds of detectors. 

The paper is organised as follows: we estimate the 
SGWB created by galactic and extragalactic exoplanets 
in Section II. Possibility of its detection is considered in 
Section III. We conclude with discussions in Section IV. 
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II. ESTIMATION OF SGWB 


B. Back of the Envelope Calculation 


A. Formalism 


The waveform and flux of GW from two point masses 
in a Keplerian orbit has been calculated very pre¬ 
cisely [49]. If a two-body system of masses M and m 
are moving in Keplerian orbits around the centre of mass 
with effective semi-major axis a, the total GW power ra¬ 
diated by the system (a.k.a. GW Luminosity), averaged 
over one period of the elliptic motion, is 

^ + m) 

^0 - 

The average frequency of this wave is /o = wq/tt, where 
Wo = V G{M + m)/a^ . (2) 


If the system has an eccentricity e then the total power 
emitted increases and becomes 
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In the eccentric case the radiation is no longer monochro¬ 
matic. The total power radiated in the n**' harmonic, at 
a frequency / = uwo/tt, is given by 
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The energy density received on earth from a source of 
luminosity L is p = L/AirdPc. We express the energy 
spectrum in terms of the usual dimensionless quantity 
^Gwif) defined as the energy density per unit loga¬ 
rithmic frequency interval in the units of the average 
critical density required for a spatially flat universe. 
Pc — 3 HqC^/SttG, where Hq := lOO/iioo km/s/Mpc is 
the Hubble constant at the current epoch. 
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The above quantity can be converted to characteristic 
GW strain [6, 29] through the equation 


hcif) = 1.5 X 10-18 y/ifooHcw)/)/"', (7) 

which is easier to compare with experimental sensitivity. 


Before indulging into the detailed numerical compu¬ 
tation of the background, we first obtain a back of the 
envelope estimate for a simple system with parameters 
reasonably close to the actual observed values presented 
in the next section. Gonsider a uniformly dense spherical 
galaxy of radius R ^ 17kpc with A ~ 2 x 10^^ stars all 
of the same mass M = IMq and each one with a planet 
of Jupiter mass m = 10“^ Mq in circular orbit with a 
period of 10^ — 10® sec. If the planets are distributed 
uniformly in logarithmic frequency interval, the average 
number density of the planets per unit frequency is given 
by A//In 10. Since the GW frequency is twice the or¬ 
bital frequency, the total GW energy density per unit 
logarithmic frequency interval received at the centre of 
this galaxy is given by 

dPGW / Ljf) A ^^^^ 2 ^ ^NLjf) 
din/ c Jq R^f in 10 JTrA^clnlO’ 


where L{f) is the average GW luminosity of one planet 
with wave frequency / (i.e., orbital frequency //2). For 
extreme mass ratio binaries such as these, where M 3> 
m, L{f) « (32G^/®/5c®)M‘*/®TO^(7r/)^°/® [see Eq. (2)]. 
Putting all together one gets 
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and in terms of characteristic GW strain 


hcU) = 1 X 10 
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If instead of a spherical galaxy, we take a disk of the same 
radius and of thickness 600 pc, the estimate doubles 

This hc{f) is overlaid on Fig 2 with a dashed line along 
with more detailed numerical estimation presented in 
the next section. Though simplistic, this calculation, 
matches the numerical estimate, including the power law 
index, reasonably well. 

Few more comments are in order here. The numbers 
used in the back of the envelope calculation are motivated 
from the distributions obtained from real data and the 
spectrum shifts by orders of magnitude if the numbers 
are not reasonably precise. So, strictly speaking, this 
is not a back the envelope calculation. The reason for 
choosing this frequency range of 10“® — 10“^ Hz for this 
calculation was that the histogram of the orbital periods 
(not shown in the paper) peaks in this range, so the esti¬ 
mation error will be minimal. Although the background 
seems low in this range, it has a rising spectrum and only 
full numerical evaluation could give a full spectrum and 
hence analyse the detectability of the background. 















3 


C. Radiation from confirmed exoplanets 

In the last few decades an extraordinary number of 
exoplanets has been discovered. The number of known 
exoplanets is growing rapidly. We collected the orbital 
data of known exoplanets from the Exoplanet Orbit 
Database [50] . Among the 1499 confirmed planets in the 
database, the parameters required to calculate GW emis¬ 
sions are available for 596 planets. Most of these planets 
are within 100 parsecs from earth. We assume the pa¬ 
rameters semi-major axis and eccentricity do not change 
significantly over the observing period. 

We then calculate the power radiated by each plane¬ 
tary system and the frequencies in which the power is ra¬ 
diated using Eq. (2-4). The total spectrum is calculated 
by adding the average flux received from all the plan¬ 
ets. Though GW radiation is not isotropic from a single 
source, while averaging over a large number of sources 
with no preferred orientation, it is reasonable to assume 
isotropic emission. The characteristic strain spectrum, 
obtained from eqn. (7) is shown in Fig. 2 by black circu¬ 
lar dots with errorbars. 


D. SGWB from galactic exoplanets 

The main aim of this paper is to estimate the back¬ 
ground from all the planets of the Milky Way galaxy, 
which we do via simulations. Studies of exoplanets sug¬ 
gest that there might be one or more planets per star 
in our galaxy [33]. We first compute the distributions 
of parameters of the detected exoplanets, taking into ac¬ 
count the errorbars, as shown in Fig. 1. We do not ob¬ 
serve any noticeable correlation between the parameters. 
We estimate the background by adding flux from sam¬ 
ple star-planet systems with parameters randomly drawn 
from the neighbourhood of their observed values and ran¬ 
domly placed with uniform distribution at different parts 
of the galaxy. We assume that there are 2 x 10^^ stars 
each with only 1 planet, though the actual number can be 
more. In reality, some stars have many planets and some 
have none, but since in general each planet emits at a dis¬ 
tinct frequency, it is statistically equivalent to distribute 
the planets uniformly among all the stars. Here we have 
taken the shape of the galaxy to be a circular disk of ra¬ 
dius 17 kpc and of thickness 0.6 kpc and placed the solar 
system 8.34 kpc away from the centre of the disk, where 
we estimate the background. Since the “quadrupole for¬ 
mula” used for calculating the Luminosity breaks down 
and flux appears to diverge if the distance is zero, in or¬ 
der to avoid random samples to appear too close to the 
solar system, we exclude the samples closer than lOOOpc 
from the solar system. We fill this “hole” with samples 
from actual observations, which lie mostly within lOOOpc, 
as seen in Fig. 1. The backgrounds estimated from this 
simulations are shown in Fig. 2 by red circular points. 
We discuss the detectability of this background in the 
next section. 


Mass of the Star 




Semi Major Axis 



Eccentricity 



Distance of the System 



FIG. 1: We make histograms of measured planet param¬ 
eters obtained from the exoplanet databases. The bottom 
histogram has been created from 596 observed for which the 
distances were known, out of the total 3988 planets, which 
were used for the top four histograms. Wherever eccentricity 
was not measured, it was set to zero to get a more conservative 
estimate. The simulations are done by drawing samples from 
the distributions of masses, semi-major axis and eccentricity. 
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Characteristic Strain Spectra 
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FIG. 2: Our estimated characteristic strains and sensitivities and upper limits of LISA and different PTAs are shown here. The 
black circles with error bars provide the background from the observed planets. When the observed distribution of parameters 
is used for simulating all the planets in the galaxy, we get the background from the whole galaxy denoted by solid red dots, 
which is the main result of this paper. In frequencies higher than 10“® the background is almost equal to the Inflationary 
Background [.51] for a tensor-to-scalar ratio of r = 0.05 (overlaid in solid black line). In that frequency range the signal is higher 
than that of compact binaries in our galaxy [52] plotted in solid green line). The sensitivities of different PTAs [53, 54] (plotted 
in bluish dot-dashed curves) are way above the strength of planetary GW. The 95% upper-limit line from EPTA data [30], 
corresponding to a spectral index of q = -1-2/3, is plotted with a purple solid curve. The LISA sensitivity is in black dot-dashed 
curve. The total characteristic strain for all the planets in the Andromeda galaxy is also overlaid in cyan squares. The total 
background created by all other galaxies in the universe up to z = 6 in LGDM cosmology and uniform comoving galaxy density 
is shown with yellow diamonds. 


Exoplanet detection methods have different selection 
biases, which are in turn sensitive to different combina¬ 
tions of masses and orbit size. The confirmed exoplanets 
are more likely to have low orbit size and high masses 
and this region of the parameter space is also more sig- 
nihcant for gravitational wave emission. Thus the esti¬ 
mate provided above should be close to reality. Also the 
distribution of mass of stars should be described by an 
Initial Mass Function (IMF). So we redid the simulation 
with the star mass distributed as an IMF [-55] and the 
results are almost identical to the previous case. 


E. SGWB from extragalactic exoplanets 

We also estimate the background created by the An- 
dromaeda galaxy (cyan squares in Fig. 2) ta SOOkpc 
away from us, assuming the same parameter distributions 
as the Milky Way, except for distance from the earth, 
which is practically the same for all the planets in An- 
dromaeda, and its mass, which is ^ 2 times that of the 
Milky Way galaxy. The background is ^ 1% of the Milky 
Way background strain. If one substitutes a constant 
distance r = ta in our back of the envelope calculation, 
one can show that the background characteristic strain 
is ~ \j2l% Rjrx "^1% times that from the Milky Way, 
consistent with the numerical result. Similarly, since the 
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Virgo cluster is ~ 1500 times more massive than the 
Milky Way and about ry = 16.5 Mpc aw ay, the char¬ 
acteristic strain from Virgo is ^ -^1500/6 i?/rv = 1.6% 
of the Milky Way, slightly greater than that from the 
Andromaeda galaxy. 

Finally we estimate the background for all the plan¬ 
ets in the universe by integrating over different redshifts. 
The spectrum of background from different galaxies falls 
off with distance and gets redshifted, but the number 
of galaxies increase with distance. The combined ef¬ 
fect requires numerical evaluation for different cosmo¬ 
logical models. We apply a procedure similar to that 
of Mazumder et al. [18] for the similar case of exoplan¬ 
ets. If each galaxy emitted J{f)df amount of energy 
in the frequency range f to f + df per unit time, per 
unit Milky Way Equivalent Galaxy (MWEG) mass, the 
integrated background spectra is given by 


^Gwif) 


^ r r , 8GJ{f{l + z))n{z) 
Jq ^ + z) E{z) 


where, n{z) is the comoving number density of MWEGs 
at a given redshift z, a(to) is the cosmological scale fac¬ 
tor at the present epoch to (generally scaled to 1) and 
E{z) is the Hubble parameter, which can be expressed, 
in terms of fractional matter and dark energy densi¬ 
ties r^A and Om, as E{z) = y/^A + (1 + z)^ for 
a spatially flat universe. Assuming that the comoving 
Milky Way equivalent galaxy density of the universe is 
n{z) ~ O.Ol/ijool^PO”^ [56], a constant, and integrating 
up to a cosmological redshift of z ^ 6 with a statistically 
isotropic ACDM model, we find that in this particular 
case the integrated extragalactic background characteris¬ 
tic strain [shown in Figure 2 by yellow diamonds] is lower 
than the Milky Way background at the higher frequencies 
and comparable or little higher at the lower frequencies. 
It is worth noting that the extragalactic background is 
expected to be isotropic, while the galactic background 
is expected to be mostly limited in the galactic plane, 
hence the latter should stand out in the extragalactic 
background. 


III. POSSIBILITIES FOR DETECTION 
A. Pulsar Timing Array 

Pulsar are very precise clocks. However the presence of 
a stochastic background can add “jitter” to the pulse tim¬ 
ing. Studying the variance of this jitter for a set of well- 
studied pulsars, a “Pulsar Timing Array” (PTA) [57], can 
provide an upper limit to the stochastic background. The 
correlations between these jitters from different pulsars 
not only improve the upper limit, but can also provide 
directional information for localised sources on the sky. 
PTA is particularly sensitive in the frequency range of 
10“® to 10“^ Hz. The expected sensitivities of different 
PTAs, in terms of characteristic strain, are overlaid in 


Figure 2 with dot-dashed curves of different shades of 
blue [53, 54]. 

Our predicted background is clearly not detectable 
with current International PTA (IPTA) [42] and does 
not seem plausible even with SKA-PTA. We overlay the 
EPTA limit on SGWB from van Haasteren et al. [30], 
which provided upper-limits for different values of the 
spectral index a. Since a power law of the form fits 
the planetary background reasonable well, we read off the 
95% upper limit on characteristic strain corresponding to 
a = 2/3 and overlay with a solid purple line. 


B. LISA and eLISA 

The Laser Interferometer Space Antenna(LISA) and 
its incarnation Evolved Laser Interferometer Space An¬ 
tenna (eLISA) are triplets of satellites arranged in an 
equilateral triangle with million kilometre arms using 
laser interferometry to measure GW. The sensitivity of 
these systems go as low as 10“® Hz [58] The expected 
sensitivity LISA mission is overlaid in Fig 2 with black 
dash-dot curve [58]. 

The predicted background is closer to the LISA sen¬ 
sitivity than any other detection method. Though the 
background peaks at this frequency range, it also drops 
drastically in this band. It could be due to some selec¬ 
tion bias in planet detection. Nevertheless, the number 
of planets among the Kepler Gandidates with orbital fre¬ 
quency above 10“® Hz is not small, so the statistical error 
in that frequency band is not very large and the estimate 
is conservative. 


IV. DISCUSSIONS 

In this paper we estimate the Stochastic Gravitational 
Wave Background (SGWB) spectrum produced by all the 
planetary systems in our galaxy. Even when there are 
nearly two billion such planetary systems in our galaxy 
the total background is still small. We however see a 
significant amount of background in the band of the pro¬ 
posed space based detector LISA (~ 10“® Hz). Though 
the total strain is still couple of orders of magnitude be¬ 
low LISA sensitivity, one can hope that a future space 
based detector will achieve such a sensitivity. 

There are planets detected which have orbital period of 
few days and also there are lower frequency planets with 
higher eccentricity which could have harmonics in this 
band. It is intriguing to note that the planetary back¬ 
ground almost smoothly merges with the galactic com¬ 
pact binary background, the later is far more dominant 
at > 10-4 Hz [52]. 

Generally the SGWB considered in literature are due 
to extra-galactic sources and the analysis is done assum¬ 
ing plane waveform. In case of galactic sources, however, 
one may need special care to account for spherical wave- 
fronts. The exoplanets of Milky Way lie in the galactic 
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plane, hence one may also have to account for the shape 
of the galaxy to improve detection SNR by performing a 
search for an anisotropic background [59]. 

Electro-magnetic astronomy has allowed observations 
of only a very small fraction of galactic exoplanets, much 
less than a millionth! Which certainly limits the aver¬ 
age statistical information on planets obtained form those 
observations. Accounting for the selection bias in detec¬ 
tion, for instance, would require an independent handle 
for better understanding of the population distribution 
of the planets. In summary, if the SGWB from the plan¬ 
ets is detected, we believe it will shine new light on our 
understanding of average statistical properties of galactic 
planet population. 
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